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Abstract 

Traditionally, the discussion about the geometrical interpretation of 
inertial forces is reserved for General Relativity handbooks. In these notes 
an analysis of the effect of such forces in a classical (newtonian) context 
is made, as well as a study of the relation between the precession of the 
Foucault pendulum and the holonomy on a surface, including some critical 
remarks about the common statement "precession of Foucault pendulum 
equals holonomy on the sphere" . 

1 Introduction 

These notes intend to analyze a classical problem in Physics, that of inertial 
forces, with the methods of elementary differential geometry of curves and sur- 
faces. Inertial forces appear in the study of the dynamics of non inertial reference 
systems, and the problem they pose is that of their interpretation in classical 
(Newtonian) terms. 

In Newton's theory, the inertial forces are a consequence of the action of 
absolute space on physical bodies; however, Mach and Einstein eradicated the 
privileged role of this absolute space as a theoretical element, and this fact led 
them to a reinterpretation of inertial forces (see the discussion in [Kin 77| . pgs 
10 — 11). Following them, we could say that Newtonian conceptions have a 
strong unilateral character (in Einstein's words "In the first place, it is contrary 
to the mode of thinking in science to conceive a thing (absolute space) which 
acts itself, but which cannot be acted upon. This is the reason why E. Mach 
was led to make the attempt to eliminate space as an active cause in the system 
of mechanics". See [Ein 56| pages 55 — 56) and this reflects on the fact that, for 
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inertial forces, the (third) Newtonian law of action-reaction is not satisfied. To 
summarize, the conceptual difficulty regarding inertial forces is that they are 
not real forces, with origin in some known physical interaction (this motivates 
the term "fictitious forces", sometimes used to describe them). 

The solution offered by Einstein is well known: he proposed a geometri- 
cal interpretation of accelerations (thus also of forces) expressed in his famous 
"equivalence principle" , according to which the inertial and gravitational forces 
can be absorbed into a non Euclidean structure of spacetime (see |Lan 87| §81, 
[Wei 72| Section 1.3). However, much care must be taken when doing this 
translation from inertial forces to geometrical (non euclidean) structures. For 
example, it is a common place (see |Mar 92j . |Sha 88\ ] that Newton equations 
in the non inertial reference system associated to the Earth, applied to the 
Foucault pendulum motion are equivalent to the holonomy on the sphere, a 
typical phenomenon of non Euclidean geometry. However, sometimes this fact 
is presented as something miraculous which calls for an explanation: geometry 
determines dynamics. But, how it is possible that the geometry of the surface 
on which it moves influence the plane of oscillations of the pendulum? We will 
argue that this kind of questions must be taken "with a grain of salt" : the 
holonomy on the sphere is not responsible of the pendulum precession, rather 
(as a result of some physical approximations), the plane of oscillations of the 
pendulum is a parallel vector field on the sphere, in such a way that a physical 
effect, its total precession in a day, coincides with the holonomy (we could say 
that both effects agree at first order) . But if we modify slightly the model to 
study the motion of the pendulum, by considering the more realistic case of an 
ellipsoid, this coincidence is no longer true. 

In brief, the whole plan of this paper is as follows. The physical analysis of 
the Foucault pendulum fits into a wider class of problems, called dynamics in non 
inertial systems; in Section[5]we present a resume of the general ideas concerning 
it. Next follows a reminder of the basic notions of Differential Geometry that will 
be used and, finally, we will discuss how to apply these results to the Foucault 
pendulum. 

Let us mention that there exist in the literature studies with a similar spirit. 
Specifically, the excellent paper of J. Oprea (see |Opr 95| ) has been a stimulus for 
this work. However, the development that can be found there can be somewhat 
confusing from a physical point of view, because it makes an explicit use of 
the inertial "forces" and, by contrast, it does not account for real forces as the 
string tension, so it can be hardly seen as a correct physical analysis. Indeed, 
this influences the crucial idea that the field of pendulum oscillations is parallel, 
as if the tension is taken into account it is no longer possible to directly conclude 
that "on the plane of oscillations of the pendulum there are no tangential forces" 
(cfr. |Opr 95| , pg 521). 

We will make use only of the machinery of the basic Differential Geometry of 
curves and surfaces, in particular, no manifold theory or Riemannian geometry 
will be needed. Suitable references are DoC 77| . [Gray 98| or |Dub 98| . 
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2 Non inertial Newtonian dynamics 



A more or less standard statement of Newton's Laws could be: 

LAW I: Free (not interacting) particles follow rectilinear uniform 
motion (including the limit case of being at rest). 

LAW II: The acceleration of a body is inversely proportional to 
its mass, and directly proportional to the external force which acts 
on it: _^ 

F = m ■ ~a . (1) 

LAW III: If a body A acts upon another body B with a certain 
force, the body B will act upon A with an equal in magnitude but 
opposed force along the line joining the bodies. 

Newton himself states that his laws are applicable just in inertial frames. A 
frame is basically a collection of physical and geometrical objects which allows 
us to assign a set of numerical coordinates to any event (Einstein described 
them intuitively as a coordinate system or reference and a clock in each point 
of space). Newton considered that such a frame with an absolute character 
exists (a very interesting discussion is given in |Rin 77| . pgs. 2 and 14-15) and 
defined the inertial frames as those moving with a uniform rectilinear velocity 
with respect to the absolute frame, an idea already present in Galileo. 

From an operational point of view, the inertial frames can be distinguished 
with the aid of the first law, but what happens with non inertial frames (NIF) as 
in the case of two observers with uniform rotational relative motion?. In these 
systems the preceding laws are no longer valid. The classical approach tries to 
continue making use of the inertial frames machinery and forces the validity of 
the second law, at the cost of introducing the so called inertial forces. Thus, in 
a non inertial frame the second law is not JQ) but 

LAW IF: The acceleration of a body is inversely proportional to 
its mass, and directly proportional to the total force which acts on 
it, including external (real) and inertial forces, F^xt and Fin'. 

Fin ^m-~a ~ Fext- (2) 

Example 1: In the particular case of a frame with uniform rotational motion 
with respect to another fixed one (this one considered as inertial). Fin takes 
the form of the Coriolis "force" : 

—To A ("tj A T*) — 27j A iTrot, 

where Tj is the rotation velocity of the rotating system with respect to the fixed 
one, r is the vector joining the origins of both systems and v rot is the velocity 
with respect to the rotating system. 

The appearance of the inertial "forces" can be seen as the consequence of 
making a change of coordinates, from a set of Euclidean coordinates (in the 
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example those of the fixed-inertial system) for which the Christoffel symbols 
equal zero, to another set of arbitrary ones without this property (see |Kop 92| , 
pg. 113). These terms involving the Tj^ are precisely those generating 
On the other hand, the Christoffel symbols are intimately related to intrinsical 
magnitudes such as curvature, so it seems plausible to ask if the effect of inertial 
forces can be understood in terms of these intrinsical magnitudes. Intuitively, 
we could say that the idea is to replace equation ((2Jl by 

F in — ^^geometrical effecf\ 

As we have remarked, for the Foucault pendulum it has been suggested that 
the geometrical effect is the holonomy on the sphere. However, as we will see, 
this must be carefully analyzed. 



3 Geodesic curvature and holonomy 

It is well known that in the case of curves in the Euclidean space the Frenet 
frame gives an M"^ basis (a reference) which is adapted to the geometry of the 
curve. Now, we will work with a curve contained in a regular surface and so we 
will consider a similar object, that is, an adapted basis not only to the curve 
but also to the surface. 

Let /S C be an oriented regular surface and c : I ^ S with c : t ^ c{t) a 
curve on it. We will call an intrinsic Frenet frame a reference along c, denoted 
i?2, £-3}, and defined as follows (the point in c denotes derivation with 
respect to t, and iV is the unitary vector in M'^ normal to the surface. The 
symbol A denotes the cross product in 



E2{t) = N {c{t)) h Ei{t) 
Es{t) ^ N {c{t)) . 

By construction, {£'i,i?2} is an orthonormal basis of T^(^f-^S, the tangent 
plane to S at c{t). If we represent by (•, •) the Euclidean scalar product in M.^ 
and by -^(i) the covariant derivative of the vector field V{t) on S (that is, the 
projection of V{t) on Tc{t)S with respect to the Euclidean scalar product of the 
ambient space), then it is clear that 

(E,{t),^{t)^ = [E,{t),E,{t))^0, 

so for each t £ I, {t) must be colinear with E2 (t) ; thus 

DEi , , I , ^ DEi , A 
-^it)^{E2it),-^it)\E2it). 
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The geodesic curvature of c at t is the proportionahty factor between (t) 
and E2{t). More precisely, let C be an oriented regular surface and 
c : / — > 5 a curve on it. Its geodesic curvature is the function 

Geometrically, Kg expresses the angle between the acceleration of the curve and 
the normal to the surface. When this angle is 0, the acceleration is tangent to 
the surface and the curve is said to be a geodesic. It turns out that, for the 
sphere, the geodesies are given by the intersection with planes passing through 
the origin, that is, the great circles. 

Let W{t) be a vector field on S parallel along c, that is a vector field such 
that -^^(0 ~ for all t <^ I. Let us call (j) the angle it defines, at each point, 
with the tangent to the curve c: 

Under these assumptions, we have a basic result (cfr. |Do(' 77j . pg. 252 or 
|One 97j . pg 338) telling us that 

that is, the variation of the angle is a characteristic of the curve c. This 
expression, is particularly suited to the practical computation of the parallel 
transport. Let ti,t2 S / be two instants, then the parallel transport along c 
between ti and ^2 mapping, denoted t(c)*J : T^^^i^-^S T^f^^^-jS, carries v £ 
Tc{ti)S to 

t{c)11{v) = \v\{Ei{t2)cOS(f>{t2) + E2{t2)sm(f>{t2)), 

where ^ 

Ht2)=Htl)- f Kg{t)\c{t)\dt. 

Jtl 

According to this formula, if we have a closed curve on S (parametrized in [0, 27r] 
without loss of generality) with c(0) = c(27r) — p, when parallel transporting a 
tangent vector v in p, making a complete loop around c, we will not recover v 
but a vector defining a certain angle with it; this phenomenon is known as the 
holonomy. 

More precisely, with the preceding conditions the angle of holonomy of c in 
p is defined as the angle between v and t{c)q'^{v), 

Example 2: Let us compute the holonomy of a parallel on the sphere 
(of unit radius). Consider the geographical parametrization, given by 

X : [-f,f] X [0,27r] ^ 5*2 

(u, v) 1-^ X{u, v) = (cos u cos w, cos u sin?;, sinit), 
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and the parallel of latitude ^o, c{t) = (cos cost, cos 6*0 sin t, sin 6*0); it is easy 
to see that its geodesic curvature is constant with value (note that the sign 
depends on the orientation of S'^ chosen) 



while 



so 



Kg{t) = — tan 00, 
\c{t)\ = cos 6*0, 



n= Kg{t)\c{t)\dt= / sinMi = 27rsin6»o. (4) 

Jo Jo 

Example 3: A model for Earth's surface more realistic than a sphere is 
an ellipsoid of revolution Ea^b (with distinct semiaxes a and b). It can also be 
parametrized geographically: 

X : [-f,f] X [0,27r] ^K,6 

{u,v) H- > X{u,v) = (a cos « cos t;, a cos u sin t;, 6 sin u), 

and a computation as in Example 2 gives for the parallel on the ellipsoid c{t) = 
(a cos ^0 cos t, a cos sin t, b sin ^o) the results 

tan 6*0 

Kg{t) = = 

V sin^ 9o + 6^ cos^ 9o 

and 

. 27rasin0o .r^ 

h= (5) 

V a? sin^ 60 + 6^ cos^ 
Of course, when a = b these coincide with the corresponding to the sphere. 

4 Parallel vector fields on the sphere and the 
ellipsoid 

Consider a curve c : / — > S' C M"^ and X : J7 C ^ 5 a paramctrization of the 
surface S, so c{t) = X{u{t), v{t)). Let F be a vector field along c such that in 
the reference {Xu{t), X^{t)} of Tc[t)S has coordinates {y^{t), y^(t)), that is, 

V{t) = V\t)Xu{t) + V^{t)X,{t) e T,(t)S. (6) 

For a vector field with that expression, it is known that the covariant derivative 
along c is given by 



DV 
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In particular, if wc consider 8 = 8"^ (the unit sphere), X the geographical 
parametrization of Example 2 and c{t) the parallel of latitude 9o, we get 



u = 9o, u = 
V = t V = 1. 



Also, for the Christoffel symbols: 



so 



DV 
~dt' 



r}^ = 0, = 0, r22=sinucosu 

r?i = 0, Tl^ = -ta.nu, T^^ = 0, 



(t) = ( V-i(t)+ sin 6lo cos 6*0^2 (v^{t) - t&neoV^{t)^ X^. 



According to that expression, a vector field along a parallel of latitude 6*0 
on the sphere is parallel if and only if its components with respect to the basis 
{Xy_(t), Xy(t)} (recall that X is the geographic parametrization) verify 

V\t)+sm0ocoseoV^{t)=O 

(7) 

If we write (3 = sm6o,x = ,y = cos^o^^ this system is equivalent to 

l + y = o 

y-l3x = 0, 
which has solution 

y{t) = B sin pt + A cos pt 

x{t) = -A sin pt + B cos Pt. 

Thus, the more general expression for a parallel vector field along the curve 
c{t) = (cos 6*0 cos <, cos 00 sin sin 6*0) (in the geographical parametrization) on 
the sphere is 

Vit) = {-Asm pt + B cos pt)Xu{t) (8) 

H 3— (B sin Pt + A cos pt) X^ (t) . 

cosyo 

Given the Christoffel symbols for the ellipsoid Eay. 



sin u cos u 



■pi 2 sin cos — a^) -pi q pi 

11 sin-^ u+hi^ cos^ u ' -^^ ' -^^ sin'' u+b^ cos^ u 

r?i = o, rf2 = -tanu, Tl2 = o, 

a similar analysis shows that the more general expression for a parallel vector 
field along the curve c{t) = (cos cost, cos sin sin ^o) (in the geographical 
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parametrization) on the ellipsoid is given by a vector V{t) = V^{t)Xu{t) + 
V'^(t)X^{t) whose components satisfy a system of equations similar to 

\ (9) 

[ cos 6*0^^(0 - sin 6*0^1 (^) = 0. 

5 Pendula, Earth rotation and holonomy 

To carry on our analysis, we begin with two basic assumptions. The first, that 
the Earth surface is spherical (and, as a consequence, the gravitational field it 
creates). The second, that the Cartesian frame with origin the center of the 
Earth sphere is (approximately) inertial. In order to do explicit computations 
we will consider the geographical parametrization of the sphere, and we will 
choose a non inertial frame associated to the parallel u = Oq as the intrinsic 
Frenet frame, that is: at each point c{t — v) the "X axis" is the direction of 
Ei{t), the "F axis" is the direction of E2{t) and the "Z axis" is the direction 
oiEsit). 



i 




Figure 1: The Frenet frame. 



Imagine the daily motion of the Foucault pendulum as the result of displac- 
ing the pendulum along a parallel on an Earth which stand still. Then, this 
non inertial frame will be rotating jointly with the pendulum. For the sake of 
simplicity, let us suppose that the mass of the bob is m = 1 and that we set 
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the time scale in such a way that 1 day equals 27r units of time, that is, the 
module of the angular velocity of the pendulum is a; = 1. If r is the vector 
giving the position of the origin of the rotating system in the Cartesian one and 
V f, a f, V rot, a rot are the velocity and acceleration (respectively) with respect 
to the fixed Cartesian axes and the rotating system then, for an observer in this 
rotating system the effective force on the pendulum is (see | Cho 95| . Sec. 10.4 
for the following development) 



where T is the string tension and ~g is the gravitational force. 

Remark: Note that the term — cu A ( w A r ) can be included in g , as g 
is only defined through the measures we make. Incidentally, this is the cause of 
the fact that the direction of ~g actually differs from the local vertical. However, 
we will assume that g is directed along Ej,{t). 

Here comes our third assumption: if the displacements of the bob have little 
amplitude, so the approximation sinC ~ C is valid (where C is the angle that 
the string determines with the vertical), and I is the string longitude, we have 



Also, by the symmetry of the gravitational field (we are supposing the Earth 
spherical). 



Fef = CLrot = Of — ~UJ A CuJ A V) — 2~iJ A Vrot- 



Thus, the equation of motion becomes 



a-rot ^ g + T - 2lU A Vrot, 



(10) 



Tx — —Tx/l 
Ty ~ -Ty/l 
T,~T 




Then, computing w A r we get 



[To A ~Vrot) 



X 



-ysin 6*0 



< {~UJ Al^rot) 



y 



X sin 6*0 



^ [To Alt rot) ^ 



-y COS 6*0, 



and the equations we are interested in are: 




Ty/l — 2x sin 00- 



Tx/l + 2y sin 6*0 
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Figure 2: Forces on the pendulum. 

For displacements of little amplitude, the tension and the weight are approx- 
imately equal T ~ 5. If we set 




/? = sin6'o, 

we arrive at the (approximate) equations 

( X + o?x = 2f3y 
\ y + a'^y = -2l3x, 

a system of coupled second order ordinary differential equations that are the 
mathematical model for the Foucault pendulum within the approximation of 
little amplitude displacements. It can be solved adding to the first equation the 
second multiplied by i, so (writing q = x + iy) we get 

q + 2ij3q + a^q = 0, 

with solution 

q{t)c^exp{-i|3t){Aexp{-t^/p^^^^)+Bexp{-t^/^P^^^)), 

a being the pendulum's pulse of oscillation. Of course, this is much greater 
than the rotational velocity of the Earth, so a ^ (3 and 

q{t) = exp{—ipt){Aexp{iat) + B exp{—iat)), 
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or, if we denote by qo{t) the solution corresponding to the absence of rotation 
(a static pendulum), 

q{t) ~ qo{t)exp{-if3t). 



Equating real and imaginary parts: 

(11) 



^(t) ] _ f cos(/3t) sin(/3t) \ / xo{t) 
y{t) ) ^[ -sin(/3i) cos(/3i) J [ yo{t) 

where (a;o(t), ?/o(0) the "static" solution for the pendulum in c(0). As is well 
known, this oscillation would always take place in a plane determined by the 
initial conditions; what is telling us is that this plane precesses when the 
pendulum is displaced along the parallel c(t), as was noted by Foucault. 

Let us check that the dynamics of the problem determines a vector field 
(associated to the pendulum oscillations) which is parallel along c, so the total 
angle accumulated after a complete loop will coincide with the holonomy of c, 
given by l@J. 

As we have already mentioned, the initial oscillation defines a unitary vector 
through the intersection of the plane in which it is contained (in M'^) and the 
tangent plane to the sphere in c(0), let us call it 
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and what we have is that {xQ{t), yo{t)) is parallel to {A, B): 

yo(i) J [b 

Now, the vector field associated to the oscillations of the pendulum (defined 
by the intersection of the oscillation plane with the tangent plane at each point) 
is 

V : [0, 27r] TS 

t ^ Vit) ^ x{t)Ei{t) + y{t)E2{t), 

where, from Hll|l . 

x{t) ^ B sin pt + A cos pt 

y{t) = ~A sin pt + B cos (it, 

so, as time goes by, the plane of oscillation actually precesses with respect to 
the plane of the initial oscillations. 

Let us note that in the geographical paramctrization the following relations 
hold: 

-^i'^'-'' ^ \\Xt(e„M - COS0O 
E2{t) - X,,{t), 
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so we can write 



V{t) = {B sin f3t + A cos (3t) —Xy{t) (12) 

cos Uq 

+ {-Asinpt + B COS l3t)Xu{t). 

Comparing with ||SJ), we see that the vector field of the oscillations is parallel 
transported, and the precession it experiments after a day (a whole loop around 
c) is 27rsin^o, according to Q). 

Remark: As a curious observation, at the equator we would not find any 
holonomy. From a geometrical point of view this is a consequence of the fact 
that the equator is the only one among the parallels of the sphere that is a 
geodesic (a great circle) and, consequently, has Kg{t) = (recall the comments 
following lO). 



6 The ellipsoidal Earth 

The purpose of this Section is to analyze what happens when we change the 
sphere model of the Earth. In order to understand it, we recall some elementary 
facts about Geodesy. 

The theoretical shape of the Earth is described by means of its normal 
equipotential surface, which coincides with the mean level sea surface (in firm 
land, the surface of a water channel) and is called the geoid. The geoid in- 
cludes the variations in gravitational potential due, among other factors, to the 
non homogeneity of the crust; this is the fundamental reason that makes the 
description of the geoid only approximate. 

The simplest approximation considers the geoid as a sphere (the ideas intro- 
duced by holonomy are based on this assumption, as we have seen), but another 
-much more precise- one is that of an ellipsoid of revolution. Indeed, conven- 
tionally, the gravimetric reductions of experimental data are carried on to a so 
called reference ellipsoid defined by the experimental values founded for some 
parameters (such as the equatorial radius of the Earth R, its total mass M or 
the flattening coefficient /) and the requirement that the ellipsoid surface be an 
equipotential surface. That is: starting from experimental data, we construct a 
mathematical model for the gravitational potential of the Earth, in such a way 
that the reference ellipsoid is an equipotential surface. The precise form of this 
potential is not of interest for us (see [Men 90j or |Tor 91| for details), but what 
we will need is to retain the fact that in the ellipsoidal model the gravitational 
field g is still normal at each point of the surface. This will be of fundamental 
importance in what follows. 

Let us return for a moment to the computations for the sphere case in Section 
^ We can see that the construction of the Frenet frame {Ei{t), E2{t), E^^t)} 
is identical in the case of the ellipsoid, as it is equation of motion UlUfl and the 
approximations following it (which are based on the fact that g is directed along 
E3{t)). Thus, for the motion of the Foucault's pendulum along a parallel c on 
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the ellipsoid we get again the equations 

V{t) = {B sin f3t + Acqs f3t)Ei{t) + {-Asml3t + Bcospt)E2{t). 

As for the parallel c{t) — X{do,t) — (acos^o cost, qcosOq s'mt, bs'mt) on the 
ellipsoid we have 



Ei{t) = 



11X4(6*0, i) II a cos 6*0 ■ 

Xu{t) 



\J a} sin^ 00 + cos^ 6*0 ' 
we have, instead of ((T^ : 

Vit) = ( B sin (3t + Acqs f3t) ^—XJt) 

a cos Oq 

+ {-AsinPt + BcosPt) ^ Xujt). 

V a2 sin^ 60 + 52 cos2 6*0 

This is to be compared with 0, which gives the expression for a parallel 
vector field along c on the ellipsoid. What we see, is that the components 

yi^^^ —A sin (3t + B COS (3t 



^2.^^ B sin f3t + Acqs pt 



\/ sin^ ^0 + cos^ 6*0 

and 

y\t) = 

a cos do 

do not satisfy lO- For instance 

V sin^ 60 + 52 cos2 Oq 

while (recall /3 = sin 6*0) 

V^{t) = ^:^:^(-ylsin/3t + Bcos/3t), 
a 

so cos9QV^{t) — sin 6*0!^^ (i) ^ 0. However, we see that the equation is true when 
a = 6, as it should be. 

To summarize: for the ellipsoidal Earth the vector field associated with the 
displacement of the Foucault's pendulum is no longer parallel, so even it does 
not make sense to speak about its holonomy. 

7 Conclusions 

We have seen that the holonomy can not be considered as the responsible for 
the Foucault pendulum precession, as changing slightly the spherical model of 
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the Earth by an ellipsoid Ea.b (which can be taken as "round" as wanted, just 
making a —>■ b) the numerical values of both effects no longer coincide (indeed, 
it is nonsense to speak about holonomy in this case), this only happens in the 
limit case a — b. From the analysis in the preceding Sections, we can see 
explicitly that the holonomy only intervenes in the precession of the pendulum 
as a first order effect: that the equation H12|) describe a parallel vector field is 
true for the sphere after making a set of physical approximations, as considering 
only little amplitudes or a spherically symmetric gravitational field. The more 
realistic situation of an ellipsoidal Earth can be seen as a perturbation of the 
spherical case, giving a first order term in which the holonomy appears, but not 
as the definitive effect. Of course, it would be very interesting to study which 
geometrical effects appear at higher orders. 
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